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Abstract 

The quantum description of Kgth propagation inside a planar waveguide is 
given , looking in particular at the behaviour of the field inside a directional 
coupler . Nonclassical effects are presented and discussed. 


Introduction 

Electromagnetic fields in optical guided wave systems are usually 
described simply by using classical Maxwell's equations, but there are cases in 
which a quantum treatment is necessary.Three purely quantum phenomena are 
known having no classical analogous; namely photon antibunching, 
sub-poissonian photon statistics, and squeezing of optical fields. If problems 
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connected with these phenomena or the evolution of photon statistics are to be 
dealt with, a quantum mechanical treatment must be used. 

Of course a inear system is not able to produce or change these 
properties but a noninear one is. For this reason, in the following the 
Hamiltonian for a nonlinear optical waveguide will be derived,and its application 
to some propagation problems will be considered.AJthough the Hamiltonian is 
quite general, emphasis is given to planar structures only and a more suitable 
approach to describe propagation fertomena is discussed. 

One of the results of having the propagation problem treated in quantum 
mechanical form is to allow for the possibility of studying how purely quantum 
effects propagate in linear systems.We will show ,for example, that a quantum 
effect as squeezing is affected by the operation of switching in a linear structure 
because of the phase changes involved in the operation. 


2. Quantization of the radiation field 

The recent experiments on nonclassica! states of light have called for a full 
quantum analysis of the electromagnetic field [1] especially in the cases of 
propagation of the fields inside dipersive media. 

We remember that the standard quantization method consists of writing the 
Hamiltonian in a given volume V , demanding periodicity in space. For 
propagating fields, the space evolution is then replaced by a time evolution, by 
finking the space and time variables by the equation z- ct. The length of the 
nonlinear medium is then replaced by an effective interaction time. Of course this 
method has two main limitations. The firt one is that , by identifying the space 
evolution with time evolution we lose one variable and this formalism can 
describe only c.w. operation : the second problem is that this procedure cannot 
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be applied rigorously to a dispersive medium, where each frequency propagates 
with different velocity. 

However we remember that by using the Hamiltonian formalsm and working in 
the Heisenberg picture , the time evolution of operator is given by 


(1) 

SO that the generator for time evolution is the Hamiltonian U. . while the 
generator for space propagation is the moment um operator fc : 

( 2 ) 


The fc? operator is related to the wave flux of the Poynting vector [2J. 

Therefore a suitable way to quantize the radiation field to describe the 
propagation fonomena is the one starting from the flux of the Poynting vector. 
This leads us to the realization that the important quantity is the flux and not , 
as usually is assumed with the Hamiltonian formalism , the energy density. 

In this way instead of quantizing the field in a large volume and demanding for 
spatial periodicity , it is necessary to assume a time periodicity T of the field , 
with the requirement that T must be forge with respect to any relevant time . 
Then instead of writing the field in term of spatial modes ( thus performing a 
Fourier analysis of the space variable z into the wave vector K m ) it is possible 
to write it in term of temporal modes ( thus performing a Fourier analysis of the 
time variable t into discrete frequencies a> m ) and space dependent 
operators. The advantage is that the temporal modes remain the same inside 
and outside the dielectric medium [3] and the space evolution of the mode 
operator can now be obtained by means of the momentum operator; moreover 
dispersion of the material can be included. 
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By the help of this formulation the expression of the electric Held is ( inside a 
dielectric ) 

E (z,t) ■ E*(2,t) + E ( z.t), 
where the cross stands for c.c. and 

fe- uo - £[ (3) 

being a (z,com) and their conjugates form a set of locaized creation and the 
annihilation operators the field frequency , i%to) the refractive index at the qd 
frequency , e o the dielectric constant and c the Igth velocity. 

The number operator for the field becomes 

N(zo.Q>m)- a + (zo.fflm) a(zo.®m) ( 4 ) 

which represents the number operator of the photons of frequency ® m 
p assi ng through the plane z ■ zo during a period T , and the commutation 
rules now become commutation at "equal space": 

fa(z;aw),a + (z\o >/>] - 5 (5) 


and the fc? operator is defined as 

fe<Z) = X(/» k/rt) a* (*®frt>a(2,®ro) ( 6 ) 

m 

where knf* n (com) ~ i« the wavevector of the field. 

c 
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3. Quantum machanical description of propagation in a planar 
wavagulda 

A planar dielectric waveguide is a medium whose dielectric permittivity 
depends on one direction, parallel to which we shall assume the x-axis (s.Fig.1). 

If this medium does not contain absorbing centers.if there is no 
amplification of radiation .and if the permittivity is weakly dependent on the field 
frequency co i ,the electromagnetic field inside the guide is expressed in terms 
of normal modes in the following form [4,5] 


A (r) - o f (x) exp( I ft .r), (6) 

where 0 j is the wave-vector with components y and z (0 j.r- k y y + k* z) of the 
j-th mode propagating inside the waveguide and f (x) is a function dependent 
only on x, defined over all space, and determined by the waveguide structure. 
Therefore each guided mode is defined by a ft vector at each coj frequency. 

From the quantum thery point of view ( as pointed out in the previous 
paragraph ,if the operator describing the field mode in a free space is given by 
( note the operator is the one that obeys at the equal space commutation rules) 

au,avi * ati)e~ ***» a*0) a~ lk *e~ (7) 

in a guided strucure it can be described as 

a(z,(oi) = atZ)6~ at 0) e~ ** (8) 
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and eq. (6) written In operatorial form becomes 


(r,0 = la(z,ciy) ®xp(“ tyO + c.c.J (9) 

i 

where A o is a constant . 
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Rfl.l A planar waveguide 

We would like to remark that in the case of dispersionless material the time 
evolution of the field operators( Heisenberg equation) is the same as the one in 
space , i.e. the Hamiltonian operator and the momentum operator 
approach provide the same results. This remark is particularly important when 
we consider the quantum treatmen of a guided mode inside a guiding 
structure, due to the fact that in general we propagate different spatial modes 
of equal frequency and we are not obliged to take into account dispersion of the 
medium if we assume a c.w. propagation. 

The same kind of consideration is still valid if we study the propagation of single 
or more modes inside a nonlinear planar waveguide with third order nonfinearity. 
These cases have been extensively studied in the paper [5]. 

In the following we analyze the case of propagation in a directional coupler , 
which is one of the most interesting guiding devices , very important from the 
point of view of its switch properties. 
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4. Behaviour of a linear directional coupler when nonclassical states are 
Involved (squeezing) 

The directional linear coupler consists in two adjacent and parallel 
waveguides [ 5,6 ](channels ).When radiation goes through the 

structure, exchange of power between the channels is possible because of the 
evanescent field which is present in the region between them. 

In the frame of classical theory the coupler is studied by using the coupled 
mode theory [5,6 ],in which a perturbation polarization responsible for the 
coupling contains the refractive index of the guides. Complete power transfer 
occurs in a distance L - ( x/2)K , where K is the coupling constant 
determined by the refractive indices of the structure; if the detuning parameter 5 
is zero, that is in the case of complete phase matching [6 ] .being 

5 = (10) 

where Pa and P& are the wavevectors of two modes of equal frequency 
propagating in channel a and b respectively .If 5 is not zero the maximum 
fraction of power that can be trasferred is proportional to 

iKF- 


From the classical equations for the complex amplitude for the directional 
coupler, in the frame of the coupled-mode theory, with obvious generalization we 
have the following Heisenberg equations for the operators, 

( da/dz - -iKbexp(i2 8z), 

/ A A 

[ db/dz ■ -iKaexp(-2i5 z), (11) 
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where a and b are the field annihilation operators in channel a and b 
respectively .For the sake of simplicity we neglect damping terms because we 
are interested in the coupling effect onlytfhis Is a good approximation in the 
region of low temperature and optical frequencies. Without taking into account 
dispersion the set of eqs. ( 1 1 ) is the same that we can write startimg from the 
momentum operator with the substitutions t -> z/c and fa -> c & 

In this way we get the following solutions of eqs.(1 1 ) 

a= C a ao+G a bo h=Cbbo+Gbao ( 12 ) 

where a> and bo are the input annihilation operators and 
C*= eP*[cos(yZ)- /54sin(yz)] 

Co= e~ ®^cos(Y 2 ) + / 5 Af sin(Y-z)] 

Ga = - i sin(Y*) e® r 

Gb= - i sin(yz) e - ®* 
where y 2 * f^ + 6 2 . 

To study the propagation of nonclassical field through the structure we use the 
following characteristic function 


Cn(P) = Tr(p exp[pa + ] P*a ] = exp[ - M| p | 2 + l* (S^p 2 + Sp* 2 )+ p lV - p W] 

(13) 


which is able to describe a field which is not a pure coherent or squeezed 
state f but has simultaneously squeezed .coherent and chaotic features [ 7 ].ln 
eq.{ 13 ) W» Wexp(i$ ) is the coherent signal, and M and S are related to the 
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noncoherent part of the fiekf.So for the vacuum state we have M- S- W« 0 
;in the pure coherent state M-S-0, and for the chaotic field W-S«0. 

The state is a pure squeezed state if 

M« Q- 0.5((4 S +1) *1) 


S- exp(i <t> ) cosh(r) sinh(r) 


(14) 


r being the squeezing parameter [13 ]. 

A mixed state is given from a superposition of a pure squeezed state with 
coherent signal W and a cahotic field described by the normally ordered 
characteristic function given by eq.( 13 ) if 
M- Q ♦ N, 


where N is the noise photon number. 


(15) 


We shall suppose that the input statistics of light in both modes can be 
described by the normally ordered characteristic function ( 13 ).Putting solutions 
(12 ) into eq.(13 ) we can see that the truncated normally ordered output 
characteri- Stic functions will have the same functional form as the input ones with 
new terms 

Mi> A4 0) I C2| + a 4 0) | Gal 2 

Sa- Si 0) C| + S& 0) Gi 

W a - H& 0) C a + M& 0) Ga (16) 

where the superscript (o) labels the input quantities, and similar expressions can 
be found for the b-mode,by interchanging the subscripts a and b. 
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We are interested in finding expressions for the variances < ( A Q ) 2 > and 
<( A P ) 2 > .(where fc}= a + a + and £*= -/( a-a + )) .In ref.(l4) it 

is shown that 

<(AQ) 2 > =1+2M+ S+ & 

<(AP) 2 > =1+2 M- S- & ( 17 ) 

Several interesting cases can be considered which depend on the way the 
coupler is feeded. 

Let us suppose first that a pure squeezed state enters channel b and a 
coherent state (or vacuum) channel a. It can be shown in this case that for 
L= * /2K and 5= 0 we have 

< (A £?) 2 >a ■ < (A h^bO, < (A P) 2 >a = < (Afc^M , 

< (A £?) 2 >a = < (A ty>b= 1. 

(18) 

This means that at the output of channel a we have an opposite squeezing than 
at the input of channel b .while the output in channel b shows no squeezing. 

A related situation is obtained when two opposite squeezed fields enter the 
two channels in the same conditions as in the previous case.ln this case 
squeezing is preserved in both channels because the field entering channel a 
comes out of channel b with opposite squeezing and the same happens with 
field entering channel b.At intermediate lengths of the coupler the squeezing is 
not completely preserved. 

A very interesting role is played by the detuning parameter . In general [14] 
If 5#o some noise is added to both channels and squeezing is reduced, and 
for some special values of 8 , noise is absent. Let us consider for example 
the case in which 5 = V3" K. . In this case y « 2k and for the same coupler 
length L='I)'/2k if a squeezed field entered channel b and a coherent or vacuum 
field channel a then we have at the output 
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( 19 ) 


< (A &) 2 > a = < (A fy*>a= 1. 

< (At?) 2 >& = < (A 

< (A P) 2 >b * < (Apfi> bO- . 

We can see that changing the detuning parameter from zero to V3K 
switches from one channel to another.This result is rather interesting for the 
purpose of measurement.The squeezed state is detected by interfering it with a 
coherent reference light and looking at fluctuations.The switching behaviour just 
described allow to preserve both the squeezed state and its reference beam. 


5 Directional and contradlrectlonal coupler with modes with small 
different frequency propagating Inside 


We have studied the problem of propagation of radiation in a coupler 
assuming two different frequencies inside the channels , with the hypothesis that 
each channel can support one only guided mode : this is possible if the two 
frequencies are quite similar. In general for a coupler the more realistic 
description of the field propagating inside ail the structure is the one which takes 
into account the superposition of the single modes propagating in each channel 
( so called supermodes); in the case of different frequencies this approach is 
particularly convenient and it is the one that we have adopted but in its 
quantum analogous, i.e. introducing this concept in the statistical dependence of 
the modes supported by the structure. 
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The motion equation of the operator describing the propagation inside the 
structure is given by ( Heisenberg form) 

^ - jj [fy«, few] , (20) 

where ] is the index corresponding to the mode ( b 1,2) and this equation is 
related to the momentum operator G , that for this case is given by 

fe( 2 ) = (z) fy*) + h (X* *i + (*) «(*) + h c ■ )• ( 21 ) 

b i 

where kp p/= is the mode wave vector , x is the coupling constant , which 

depends on the refractive index n(<aj ) distribution inside the coupler . It is very 
interesting to observe that the G operator looks like the one of a second order 
nonlinearity for a bulk material. 

Using the approach of the supermodes we can describe the two fields of 
different frequency supported by the strucure as 

1 1> = 2 _v *(| 1w>+ mi 1 loop* + U2i|2«rf> ) 

| 2> = 2 _V6 (|2/n>+ U12 | loop* + U22| 2ouP> ) (22) 

being uy,* a function related to the transformation law of the coupler , 
containing all the informations about the structure , such as the coupling 
constant , the detuning parameter , etc. ( see functions C a* and G a,b of the 
previous paragraph ). 

To follow the statistics of the field we start from the characteristic function ( the 
antinormal one) from which it is possible to derive all the factorial moments and 
the photon counting distribution. As in the previous paragraph we suppose the 
input state is a superposition of a coherent state and noise , including 
squeezing . 
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Due to to the hypothesis of the supermode we can write the charateristic 
function for all the fields Pi and fe . Ca (Pi.fe) , but we can follow also the 
behaviour of each separate mode C*(P/) [ 8 ] : 

Ca„(Pi.P 2) a [— Bjin | Pyl + g (tfb $ ^ c.c.) + Py ~ c.c.) + 

(noise) ( squeezing) ( coherent) 

+ (-Si2mPiP2+ C?2*,PiPz + c.c.)]}. 

( interference of noise ) (23) 

The output characteristic function is of the same form as the input, where all the 
features of the coupler are inside the B and C coefficients of the eq.(23). 

Several cases of inputs states have been studied, such as coherent , two-photon 
coherent , two-mode squeezed states and an factorial moments have been 
calculated [8] finding as the detuning parametr plays a very important role on the 
evolution of the fileds : it adds addidional noise if it is non zero [8]. 

It is interesting also to follow the photon counting distribution which put into 
evidence the switch properties of the structure always starting from the hypotesis 
of supermode supported by the coupler . An example is shown in Rg.2 , where 
the detuning parameter 6 is zero, the input state in the first channel is a two 
photon coherent state in the first mode , with a " small amount of 
squeezing" ,and a coherent state in the other mode . The picture shows the 
marginal photon number distribution in the channel 2 ; at a suitably distances the 
sub-Poissonian behaviour turns super-Poissonian, which characterizes the field 
in the squeezed vacuum state: this confirms the switching of figth of certain 
photon statistics from one mode to the other one. 
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Fig.2 - The marginal photon number distribution for 8 »0 and K «i for the 
channel 2 


Conclusions 

The propagation of Ugth in a linear directional coupler can be studied without 
taking into account the dispersion of the dielectric constant until c.w. field 
propagation is considered; of course dispersion must be taken into account in 
non stationary cases and when the structure of the propagating device supports 
different frequencies. 

The ability of the coupler to switch from one channel to the other by introducing a 
p hase lag allows to change the squeezing directions, until the ^ parameter 
is of suitably values; in general a detuning different from zero reduces the switch 
properties of the coupler and adds additional noise to the propagating fields. This 
effect is in turn evident also on the photon counting distribution. 
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